ASYMPTOTICS FOR A VARIANT OF THE MITTAG-LEFFLER 

FUNCTION 



STEFAN GERHOLD 



Abstract. We generalize the Mittag-Lefflcr function by attaching an expo- 
nent to its Taylor coefficients. The main result is an asymptotic formula valid 
in sectors of the complex plane, which extends work by Le Roy [Bull, des sci- 
ences math. 24, 1900] and Evgrafov [Asimptoticheskie otsenki i tselye funktsii, 
1979]. It is established by Plana's summation formula in conjunction with 
the saddle point method. As an application, we (re-)prove a non-holonomicity 
result about powers of the factorial sequence. 



1. Introduction 

For a,b,a > 0, the series 

oo „ 



(i) K7(*) 



4< T(an + b) a 



n=0 



defines an entire function of z. The Bessel function Io(2y / i) = X)^Lo z "/ n ' 2 an< ^ 
the (generalized) Mittag-Leffler function E a ${z) = Y^=v z n /T{an + b) are special 
cases. The function ([TJ has the order 1/aa; see, e.g., Titchmarsh [131 Example 8.4], 
for a = b = 1. His argument trivially extends to a, b > 0. 

Le Roy [8] has obtained the asymptotics of F£Y ( z ) as z tends to infinity along 
the real line. This result (a special case of formula ^ below) appears also in 
Hardy [6, p. 55]. Le Roy establishes it by converting the sum into an integral and 
then applying the Laplace method. Note that the latter can also be applied in a 
direct way: For real z, the summands of 

( 2 ) *i ( ?(') = £5= 

n=0 

are positive and concentrated near n ~ z 1 ^", and the Laplace method is easily 
carried out. The contribution of the present note is an extension of the asymptotics 
to a, b > and complex values of z, which is presented in Section[2] (For a = b = 1 
and complex z, this question is also discussed in Evgrafov's book [2]; see the end 
of Section [2J for detailed comments.) 

As a small application, we prove in Section [3] that the (possibly formal) series J2J 
is not D-finite [T2J for all a € K \ Q. (This is a special case of a known result pQ.) 

We note in passing that the function (J2j) satisfies the integral relation 



(3) / 
Jo 



- t '*Flf l \t)dt = zF[°?(z), z^Q. 
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Indeed: 



Jo n]a+1 £S " !Q+1 ^ 



00 _ /•oo 00 +1 



n=0 u n=0 

2. Main Result 

Theorem 1. Let a, a, b > and s > be arbitrary. Then, for z-iooin the sector 

!-^acrK — e < aa < 2 

(2-iaa)7r-£ 2<aa<4 
4 < aa, 

we Ziawe the asymptotics 

(4) «w4(2^^e^, 

fly a 

Applying the Laplace method directly does not work for non-real z; the absolute 
values of the summands in (fl) are peaked near n ps a -1 ^ 1 / 00 , but it seems that 
one cannot balance the local expansion and the tails. This is caused by oscillations 
in the summands, which can be dealt with by shifting the problem to the asymp- 
totic evaluation of integral. The Laplace method then succeeds, after moving the 
integration contour through a saddle point located approximately at a~ 1 z 1 / aa . 

Lemma 2. Let a, a, 6 > and e > be arbitrary. Then, as z — > 00 in the sector 
I ar g( z )| _• max{0, (2 — ^aa)Tr}, we have 

^ F S , W = rr(^ d * + 0W - 



Proof. First, let us fix a z in this sector. Put f(t) = z t /T(at + b) a . Plana's 
summation formula [7] Theorem 4.9 c] yields 

(6) e m = r mdt + 1/(0) + i r n % n ; it] dt. 

n=0 J J e 1 

To check the validity of (JSJ), we have to verify that 

lim |/(a;±iy)|e- 2 ^=0, 

uniformly for x in finite intervals in [0, 00 [, and that 

\f(x±iy)\e- 2 ^dy 



exists for x > and tends to zero for x — > 00. To do so, first note that 

= | 2 |^ e Taarg(z) < | z | x & y\ arg(z) | 

Furthermore, by Stirling's formula, we have 

|r(x±iy)""| < e ^+S)x x - a x e ( a ^/2+8)y 

for large \x + iy|, where 5 > is arbitrary. Hence 

|r(aa; + b ± ia V y a \ < e {aa+s "> x {ax + ^"K'lel™/^*, 

and thus 

p (aa+8)x\ \x 

(7) \f(x ± iy)\ < -± ' ' e<l ^ g (z)\+ aa ,/2+S)y 

w lJy yn ~ (ax + b) a i ax + b ) 
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which implies both required conditions. Finally, putting x — 0, we see from ([7]) that 
the second integral in © is O(z) as z — > oo. Since /(0) = 0(1), we are done. □ 

Proof of Theorem [7J We apply the saddle point method to the integral in ([5]) . To 
locate the saddle point, we equate the derivative of the logarithm of the integrand 
to zero, which leads to the equation 

T'(at + b) , . , . aa 

= logz - aa-^—^ = log, - aalog(at + b) + + 0(r ). 

By bootstrapping, we find that there is an approximate saddle point at 

f _-l „l/oa , 1-26 

to — a z + 

We change the integration contour to a line C that begins at and passes through tp. 
Note that this change of contour is valid for large \z\, by Stirling's formula, as long 
as | arg(to) | is bounded away from n/2. But this follows from our assumption 
on arg(z). 

The dominant contribution of the integral arises from the range 

|* — * | < l*o 

around the saddle point, where (3 is an arbitrary member of the interval ]g, |[. We 
write 

t = ta(y + l), -l<y<oo, 
and divide the integral as follows: 

z *o(»+l) 



dt = *0 / -F7-T7 — TTT— TT77 d 2/ 



c T(at + b)<* u 7-i T(at Q (y + l) + b) 



= to I / + / + / ^ 77. , T\ , ^ n dy 



-1 

(8) =:/i+I 2 + /3. 



Itol' 3 " 1 J\t \e-i) r («*o(2/ + l) + 6) a 



First we investigate the central integral /a- From Stirling's formula we get the 
following local expansion of the logarithm of the integrand: 

(9) (t (y +1)) log z~a\o g r((at (y+l) + b) = 

to log z — aato log to + cux(l — log a)to + a{ \ — b) log to 

+ a(h — b) log a — a log V2tt — ^aat y 2 + o(l). 



Since 



I* II 3 - 1 It 1)3-1/2 

exp(—^aatoy 2 )dy~ — / exp(— ^aato\to\~ 1 u 2 )du 



\to\f 3 - 1 ' ' ' ' Vl*o| "'-Itol' 3 - 1 / 2 

1 



exp(— i aato | to I 1 u 2 )du 



2tt 

aato ' 



the central part I2 thus satisfies 

(10) h ~ to2 to to aato e aQ(1 " losa)to (at ) a( ^ b) (27r)- a ' /2 y ,r ^" 



aato 



1 l-a Q-26Q + 1 1/a 

— =(27r) 2 z e 



aya 

This is the right hand side of 
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It remains to show that the integrals I\ and I3 are negligible. By our assumption 
on arg(z), there is a constant c\ > (independent of z) such that 

|S*(*o)| < cxK(to), 

hence 

l*o| 



(11) 



K(*o) > 



Now divide the integral I3 further into 




C2 *0 



^*o(y+i) 



r(ai (y + !) + &)« 



dy 



where 

u := exp(iaac2|io| 2 )■ 
By Stirling's formula, there is a positive constant cq such that 



(12) 



T(at + b) a 



< cq exp 5R( (at + 6) log z — a(at + b) log(ai + 6) 



+ a{at + b) + ±alog(at + b)j, t G £. 



An elementary calculation shows that the right hand side of (fT2")l decreases w.r.t. \y\ 
for large |z| and \y\ > |t | /3_1 - Therefore, we can estimate J31 by inserting y = 
|io| /3_1 into ([12)) and multiplying by the length of the integration path, which is 
u — l^ol^ 1 < u - Using (]9]), and writing A = A(z) for the factor in front of the 
square root in PU|) . we obtain 

vt 

< c 



r(ai + o) Q 



y=l*ol' 3 - 1 



^ e -aato|t | 2f! " 2 /2 



- colAle-^^)!* ! 2 ""/ 2 < c \A\e- aac ^of/ 2 . 
The latter inequality follows from (fTTj) . Hence 

II31I < c u|A| • e-^l*"!^" 1 ^ = colAle" ^ 2 !* ! 2 ^ 1 / 6 . 
Now we compare this estimate with (flQ)) . Since 



the integral J31 is indeed negligible. As for / 32 , it easily follows from Stirling's 
formula that 



T(at + b) a 
for large \z\. We thus obtain 



V>u, 



/32I <t 



My 



= i c- ll </ 2 . 

Finally, the integral I\ in © can be estimated analogously to /31. 



□ 



A full asymptotic expansion can be obtained easily by pushing the local expan- 
sion around the saddle point further. 

Evgrafov [2J § 4.2] offers a similar asymptotic treatment of F\ x (z) in sectors 
of the complex plane. For a < 2 and | arg(z)| < ^air — e, his result agrees with 
ours. (Except that a factor a -1 / 2 , or p 1 / 2 in Evgrafov's notation, is missing from 
the formula.) However, he gives few details on how to carry out the saddle point 
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analysis, in particular, on how to do the tail estimates. For a > 2, Evgrafov [5] 
p. 294] appears to go beyond our Theorem[T] in that he claims (HJ (with a = b = 1) 
for any sector that stays away from the negative real axis. There seems to be a 
serious gap in the proof, though. 

To be specific, we switch to Evgrafov's notation. On p. 292, he writes that 
E^Lo ^™ n ' _1 ^ P satisfies the assumptions of Theorem 4.2.2 for all p > 0. (The 
text says Theorem 3.2.2 instead, but this is certainly a typo.) This means that 
fi(z) = T(z + l)~ 1//p should satisfy \fj,(x + iy)| < M,4exp(— Ax), for arbitrary A 
and some other constant Ma, and for x + iy in a domain D containing a horizontal 
strip that contains the positive real line. But due to the exponential decrease of 
the Gamma function towards ±ioo, this can hold only if the elements of D have 
bounded imaginary part. Then also the contours Cj^ and C~\ on p. 292 must 
have bounded imaginary parts. On p. 293, the saddle point method is applied to 
the second integral on p. 292 (over the contour C_ x ), where the location of the 
saddle point is t p exp(27rip). The imaginary part of this point is not bounded for 
large t. Hence the validity of the necessary change of integration contour is not 
proven. These remarks seem to justify another study of the problem, provided by 
the present note. 

We close the section by a possible question for future research. For complex a 
and fixed z, one might ask whether the function defined by ([5]) has an analytic 
continuation for 3t(a) < 0. (The relation ([3]) does not seem to be useful in this 
respect.) 

3. An Application: Non-Holonomicity 

An analytic function (or formal power series) is called £>-finite [15], or holonomic, 
if it satisfies a linear ODE with polynomial coefficients. An equivalent condition is 
that its power series coefficients satisfy a linear recurrence with polynomial coeffi- 
cients. 

There has been some interest recently in showing that certain series (resp. se- 
quences) are not holonomic [TJ [3] 21 [TO] E] • Lipshitz 0| Example 3.4(i)] mentions 
(without proof) that the sequence (n! a ), which satisfies 

(n + l)! Q = (n+l) a n\ a , 

is holonomic if and only if a is an integer. See [TJ Theorem 4.1] for a proof that 
there is indeed no recurrence with polynomial coefficients for a G C \ Z. 

From Theorem[TJ we can conclude the weaker result that (z), and thus (n! Q ), 
is not holonomic for a € R \ Q. Indeed, since multiplication by the holonomic 
sequence n!~L Q J preserves holonomicity, we may assume that < a < 1. Hence 
Theorem[T]yields the asymptotics of F a (z). But a function that features the element 
exp(az 1 / Q ) in its asymptotic expansion at infinity, in a sector of positive opening 
angle, can be holonomic only for rational a. This follows from a classical result on 
the asymptotic behavior of solutions of linear ODEs. (For details on this method of 
showing non-holonomicity, see [3] and [?]•) Including the parameters a and b, one 
can conclude other non-holonomicity results from Theorem[lJ but they are going to 
be weaker (in terms of parameter ranges) than the corresponding results deduced 
from the method of [T]. 

Acknowledgment. I thank Julia Eisenberg for her linguistic help concerning 
the relevant passages of Evgrafov's book [2]. 
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